function g = sumVarianceGradient(X)

[N D] = size(X);
C = createCentringMatrix(N);

g = zeros(size(X));
for d = 1:D
    g(:, d) = quadraticGradient(C, X(:, d));
end
g = g/N;
